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Abstract. Pseudo-Newtonian gravitational potential introduced in spherically 
symmetric black-hole spacetimes with a repulsive cosmological constant is tested 
for equilibrium toroidal configurations of barotropic perfect fluid orbiting the black 
holes. Shapes and potential depths are determined for the marginally stable 
barotropic tori with uniform distribution of specific angular momentum, using 
both the pseudo-Newtonian and fully relativistic approach. For the adiabatic 
(isocntropic) perfect fluid, temperature profiles, mass-density and pressure profiles 
and total masses of pseudo-Newtonian and relativistic tori are compared providing 
important information on the relevance of the test-disc approximation in both 
the approaches. It is shown that the pseudo-Newtonian approach can be precise 
enough and useful for the modelling of accretion discs in the Schwarzschild-de 
Sitter spacetimes with the cosmological parameter y = AM 2 /3 < 10 -6 . For 
astrophysically relevant black holes with y < 10 -25 , this statement is tested and 
shown to be precise in few percent for both accretion and excretion tori and for 
the marginally bound, i.e., maximally extended tori allowing simultaneous inflow 
to the black hole and outflow to the outer space. 



PACS numbers: 04.20.-q, 04.25.-g, 04.70.-s, 95.30.Lz, 98.62.Mw, 98.80.Es 

1. Introduction 

Cosmological observations of distant la-type supernova explosions indicate an acceler- 
ating universe. Starting at the cosmological redshift z ~ 1, the accelerated expansion 
should be generated by some appropriate form of the so-called dark energy [TJ [2]. 
These results are in accord with a large variety of cosmological tests including grav- 
itational lensing, galaxy number counts, etc. [3J. The recent detailed studies of the 
cosmic microwave background (CMB) anisotropics indicate that the energy content 
of the dark energy represents ~ 74.5% of the energy content in the observable uni- 
verse, and the sum of energy densities is very close to the critical energy density p cr it, 
corresponding to almost flat universe [H [5] . 

A large variety of possible candidates for the dark energy is discussed in these 
days. First of all, there is the standard possibility represented by the cosmological 
constant A. Its Lorentz invariant form enables interpretation in terms of a ground 
state or vacuum energy of quantum fields [6l [7]. The energy density p\, which 
can be associated with the cosmological constant, remains unchanged during the 
cosmic expansion, and its pressure to energy density ratio (equation of state) is 
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w = Pa/pa = ~ 1- Further, there is a variety of scalar fields evolving outside of their 
energy minimum, called quintessence, which possess a time varying energy density 
and equation of state with w < —1/3 [5J. Such a scenario can be realised by light 
scalar field coming from modified f{R) gravity [5], string- inspired cosmologies [10] , 
cosmology with extra dimensions [IT] , or by fc-essence being a scalar field with a non- 
canonical kinetic term [12j [13] . Similar behaviour is exhibited by the coupled dark 
energy, i.e., a scalar field coupled to the dark matter. For example, Chaplygin gas or 
its generalization called quartessence explain both dark energy and dark matter from 
an unified physical origin [14) . The holographic dark energy treats vacuum energy 
limited by a holographic conception on degrees of freedom [15] . Modified-gravity 
origin of the cosmic acceleration comes also from multidimensional theories of the 
braneworld-type cosmology [16 . All of these models have the equation of state with 
the parameter w varying during the cosmological expansion. Another possibility is the 
evolving vacuum energy density (cosmological constant) introduced in the framework 
of superstring quantum gravity [17] . On the other hand, the present accelerated state 
of the Hubble flow can be explained in the framework of the Tolman-Bondi models 
[j~8] as a consequence of inhomogeneous distribution of matter in space. 

Recent observational data indicate that the allowed equation of state for the 
dark energy is very close to the case of a repulsive cosmological constant A > 
with pa = 0.745/9 cr it. Therefore, it is relevant to consider the influence of A > in 
cosmological and astrophysical phenomena. The cosmological relevance of A > is 
discussed in the standard textbooks [IH [20] . 

We have studied possible effects of A > in astrophysically motivated problems 
[21"] . concentrating on radial geodesies [22] and their implications in cosmological 
models [23J [2H [35] , equilibrium (static) positions of spinning test particles [55] [S7J [35] , 
and astrophysically most important disc accretion [29] . investigating circular orbits 
of test particles in the equatorial plane [30] [31] and equilibrium configurations 
of barotropic perfect fluid tori in black-hole (BH) and naked-singularity (NS) 
Schwarzschild-de Sitter (SdS) and Kerr-de Sitter (KdS) spacetimes [521 133] ■ Elliptic 
integrals of more general motion of test particles are described in [33] [35J [35] . The 
velocity profiles of the geodetical discs and perfect fluid tori in the locally non-rotating 
frames were studied in [37] [38j [39] [40] [41] [42] ■ Further, the test particle and perfect 
fluid properties were treated also in the framework of the optical reference geometry 
[43] 144] [45] allowing introduction of inertial forces in the intuitive Newtonian way 
[46] . Motion of photons in the Reissner-Nordstrom-de Sitter and more general Kerr- 
Newman-(anti-)de Sitter spacetimes is discussed in the papers [17] 25]. The optical 
phenomena were treated, e.g., in [49l [50]. [5l l [52] , [53 1 [54 ] [5% ] [56 ] [57]. 

Self-gravitating fluid solutions of the Tolman-Oppcnhcimcr-Volkoff equations with 
a non-zero cosmological constant (ATOV) were first studied in [55] for spherically 
symmetric stellar-like configurations with uniform energy-density distribution and 
generalized for cosmological solutions in [53] • Solutions of the ATOV equation for 
polytropic and adiabatic self-gravitating configurations were investigated in [60] . 
Stability of such stellar- like solutions was discussed in [61] and [62]. The self- 
gravitating spherically symmetric polytropic configurations with the dark energy in the 
form of quintessence give interesting variety of spacetime singularities [63] that enlarges 
the singularities occuring in the models based on the ATOV equation. Here we restrict 
our attention to the test fluid configurations representing physical properties of black 
hole backgrounds with the cosmological constant. To some extend our results could 
reflect also the properties of black hole solutions of the f(R) gravity (for review see, 



Pseudo-Newtonian and general relativistic tori in SdS spacetimes 



3 



e.g., [63]) where the cosmological constant term is present. Those solutions generally 
include other parameters whose influence is worth to study. 

The influence of A > on the BH-spacetime structure can be properly represented 
by the dimensionless cosmological parameter y = AAf 2 /3(3 For SdS black holes 
admitting existence of stable circular geodesies, which are necessary for the existence 
of accretion discs, the cosmological parameter y < y ms ,c = 0.000237 [3D]. Cosmological 
tests using the supernova magnitude-redshift relation and the measurements of CMB 
fluctuations [H [S] imply A ~ 10 _56 cm -2 , and thus very low present values of y for 
astrophysically realistic black holes. The most massive currently observed black hole 
seems to be harboured in the quasar TON 618, having the mass of 6.6 x 10 10 M Q 
[65l [66] . In this case, y = 4.1 X 10~ 25 . For primordial black holes in the very early 
universe, with expected high values of the effective cosmological constant, the values 
of y can be much closer to j/ mSjC . Considering the electroweak phase transition at 
Tow ~ 100 GeV, we obtain an estimate of the primordial effective cosmological constant 
A ow = 0.028 cm -2 , while at the level of the quark confinement at T qc ~ 1 GeV we 
obtain A qc ~ 2.8 x 10~ 10 cm -2 and consequently higher values of y [29l 132] . 

Geometrically thin (Keplerian) discs are characterized by the quasi-circular 
motion of test particles along stable circular orbits (geodesies) in the equatorial plane 
of a given spacetime. The inner edge of the Keplerian disc thus corresponds to the 
innermost stable circular orbit. In the SdS and KdS backgrounds, there is also the 
outer marginally stable circular geodesic which puts a natural limit on the maximal 
extension of Keplerian discs in the spacetimes with A > [29 . 

Equilibrium configurations of barotropic perfect fluid tori are described by 
equipotential surfaces of the 'gravito- centrifugal' potential W, which correspond to 
the surfaces of constant pressure in the fluid. Analyzing properties of the equipotential 
surfaces in the case of marginally stable tori, i.e. tori with uniform distribution of the 
specific angular momentum £(r, 9) = const, orbiting in the SdS and KdS backgrounds, 
three astrophysically relevant phenomena were found [32J, [33] : (i) the outer cusp in 
the structure of closed equipotential surfaces, corresponding to the outermost edge 
of the torus, (ii) strong collimation of open equipotential surfaces near the axis of 
rotation of the disc that could have some relevance for the collimation of huge jets 
outside some giant active galaxies, (iii) for current value of the effective cosmological 
constant Ao ~ 10~ 56 cm -2 , there is a dimensional coincidence between extension 
of large galaxies and maximal extension of perfect fluid tori orbiting supermassive 
black holes with mass in the range of (10 6 -10 10 ) M Q . For this interval of masses and 
current cosmological constant Ao, the cosmological parameter y is in the range (10 -34 - 
10 -26 ). Important study of stability of fluid tori in the SdS spacetimes against the 
so-called runaway instability was realized in [67] , suggesting a strong stabilizing effect 
of the outflow of matter through the outer cusp of toroidal thick discs. Therefore, the 
cosmological constant strongly influences the structure of disc configurations (both the 
geometrically thin and thick) introducing quite naturally the outer edge of accretion 
discs. 

It should be stressed that all of the relevant effects of A > on the BH physics 
are well expressed in the SdS spacetime, since the effects of dragging of inertial frames 
caused by the BH rotation are concentrated in the vicinity of the BH horizon, where 
the influence of A > can be abandoned for realistic values of the BH mass and the 
present relict cosmological constant Ao. Of course, the efficiency r/ of the accretion 

| We use the geometric system of units (c = G = 1) hereafter. 
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process is influenced by the frame dragging in the innermost part of the disc, where the 
KdS spacetime structure is relevant. For near-extreme KdS black holes characterized 
by the cosmological parameter y <C 10~ 6 , r\ = 0.4, which is much higher in comparison 
with the SdS black holes where r\ = 0.059. Therefore, in studying the large scale 
properties of disc structures, investigation of the SdS spacetime is quite sufficient, 
only the accretion efficiency has to be given by the KdS spacetime structure. 

One of important problems that could be hardly solved in the framework of full 
general relativistic (GR) approach is the influence of A > on the structure of self- 
graviting discs. For the current stellar-mass and super-massive SdS black holes, the 
strong gravity near the BH horizon r^ w 2M weakens with the distance from the 
black hole growing and, at r > M, the structure of the disc can be well described 
by the Newtonian theory. However, the Newtonian theory starts to lose its validity 
behind the so-called static radius r s = y _1 / 3 M, where the repulsive effect of the 
cosmological constant starts to be relevant up to the other strong gravity region 
near the cosmological horizon r c w j/ _1 / 2 M. We expect that the pseudo-Newtonian 
(PN) approach, based on the replacing of the Newtonian gravitational potential by 
an appropriate PN gravitational potential in the framework of the Newtonian theory, 
could be succesfull in order to describe properly some GR effects. The PN approach 
enables to use directly the standard techniques developed in the framework of the 
Newtonian physics, e.g. the Newtonian discoseismology [53], and include the repulsive 
cosmological constant and some other relativistic effects into account. 

These are the reasons that led us to introduce, in close analogy with the standard 
approach of Paczyhski and Wiita [69] related to the Schwarzschild spacetime, a PN 
gravitational potential for the SdS spacetime, and test its accuracy by comparing 
its predictions with the exact GR results [70]. We focused on the properties which 
play a crucial role for a geodesic motion and accretion processes in thin accretion 
discs, e.g. locations of the marginally stable and marginally bound circular geodesies. 
We concluded that the PN gravitational potential introduced in [70] seems to be 
suitable for the thin disc description in astrophysically realistic SdS spacetimes with 
the cosmological parameter y < 10 -25 [70] , 

In the present paper, we study the PN description of perfect fluid tori orbiting the 
SdS black holes that approximate realistic thick accretion discs, and compare the PN 
results with the exact GR results. The influence of the discs on the spacetime geometry 
is assumed to be negligible. We study both the global geometrical structure of toroidal 
configurations and internal physical properties; the latter under simplified assumption 
of adiabatic tori. In section [5] we introduce the PN gravitational potential for the 
SdS spacetime and briefly summarize its significance for geodesical motion. In section 
[31 equilibrium configurations of barotropic perfect fluid tori are calculated using both 
the exact GR and approximate PN approaches. Namely, we calculate loci of cusps of 
the equilibrium configurations, shapes of the barotropic tori and depths of the related 
potential wells; the GR and PN results are compared. In section^ we study properties 
of the adiabatic tori in both the approaches, i.e., we determine mass-density, pressure 
and temperature profiles. In section [5] we calculate total mass of adiabatic tori in 
both the approaches and give some limits on the test-disc approximation used in our 
calculations. In section[6l behaviour of thermodynamic quantities across the adiabatic 
tori and the definition of the PN potential are discussed. In section[JJ some concluding 
remarks are presented. 
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2. Pseudo-Newtonian gravitational and effective potentials 

According to a general heuristic method |71J . the PN gravitational potential can be 
defined by the relation 

f I 2 

* = J -f dr, (1) 

where r represents the radial coordinate and £ c — L c /E c is the GR specific angular 
momentum, i.e., the ratio of the conserved angular momentum and energy per particle 
mass, related to circular geodesies in the equatorial (central) plane. This definition 
is based on the Newtonian relation for the gravitational potential = / l 2 /r 3 dr, 
where l c is the Newtonian angular momentum per mass of the particle moving along 
the circular orbit[§| As shown below, this definition of the PN gravitational potential 
yields a number of advantages. Note that it also implies the well-known Paczyhski- 
Wiita gravitational potential [55] 

M 

describing with high precision the test particle motion and accretion disc structure in 
the field of any Schwarzschild black hole. 

2.1. Circular geodesies in SdS spacetimes 

Geometry of the SdS spacetime in Schwarzschild coordinates (t, r, (f>) and 
geometrical units (c = G = 1) is described by the line element 

ds 2 = - (l - ™ - yr 2 ^j dt 2 (3) 
9M \ _1 

2 A„1 i J/j/i2 , • 2 , 



1 yr 2 j dr 2 + r 2 (d6 2 + sin 2 6d<j) 2 ), 

where M is the total mass of the central object (e.g. a black hole). In the following, we 
put M = 1 in order to have completely dimensionless formulae. Pseudo-singularities 
of the line element §5§ give the loci of the BH and cosmological horizons 

Hi = cos — - — , r c — —= cos — - — , (4) 

where £ = cos -1 ('i^/iy). Both horizons are determined by the condition 

r — 2 

y = Vh(r) = —j-, (5) 

and exist for < y < y CI n — 1/27, separating the spacetime into two dynamic regions 
and one static region. For y — y cr n, both the horizons coalesce at the photon circular 
orbit located at r = 3. For y > 1/27, the horizons disappear and the corresponding 
SdS spacetimes become dynamic NS spacetimes (figured]). 

The geodesic motion in the SdS spacetime is fixed to central planes. In the 
equatorial plane (9 = tt/2), the motion is determined by the GR effective potential 

m 

V 2 s (r; y,L) =(l--- yr 2 ) ( 1 + ^ V (6) 



r 



§ The PN potential given by Eq. JTJ is denned to match the Newtonian angular momentum per 
particle mass on a circular orbit with the GR specific angular momentum. We discuss this kind of 
the PN potential definition in section [6] 
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Figure 1. Characteristic radii of the SdS spacctime in dependence of the 
cosmological parameter y. Solid curve determines the BH and cosmological 
horizons, dashcd-dotted line determines the static radius, dotted line shows 
location of the photon circular orbit at r = 3, dashed curve determines marginally 
bound circular orbits and the full thick curve marginally stable circular orbits. 
Circular orbits exist in the shaded region only, its dark part corresponds to the 
stable ones. Black-hole spacetimes exist for y < 1/27 = 0.037, those admitting 
motion along stable circular orbits exist for y < 12/15 4 = 0.000237. 



where L is the angular momentum per unit rest mass of the particle. The motion is 
allowed in the region where particle's energy per unit rest mass E > V e g. Circular 
geodesies are determined by local extrema of the effective potential. At a given radius 
they are characterized by the constants of motion 

L,= {r(l-Vr 3 )]" 2 (l- 3 -y" 2 ■ (8) 
Corresponding GR specific angular momentum then takes the form 

l c — f „ SN • yd) 
(r — I — yr a ) 

Circular orbits of test particles exist in the region 

3 < r < r s = JT 1 / 3 , (10) 

where the lower limit is given by the photon circular orbit, which is located at r = 3 
independently of the cosmological parameter y, see (3D] for more details, and the upper 
limit corresponds to the so-called static radius, which is the place where a geodesic 
observer with just time-component of its 4-momentum being non-zero can exist. At 
the static radius, the gravitational attraction of the black hole is just balanced by 
the cosmic repulsion represented by the cosmological constant, i.e., test particles 
with zero angular momentum feel no 'force' there, similarly to the radial infinity 
in asymptotically flat spacetimes. 




Figure 2. PN gravitational potential tfi given by relation \12\ for three values of 
the cosmological parameter y: 10 — 4 (solid), 10~ 5 (dashed), and (dotted). 



2.2. Pseudo-Newtonian gravitational potential 

Constructing the PN gravitational potential for the SdS spacetime, we have to reflect 
appropriately both the gravitational attraction of the black hole and the repulsive 
effects of the cosmological constant. The PN gravitational potential is defined by 
relation {!]) in which the specific angular momentum of a particle on circular orbit is 
given by relation After integration, we obtain a general form of the potential t/j, 
i.e. 

^ = -7T, ^ JT+C, (11) 

2(r — 2 — yr 6 ) 

where C is an integration constant having no physical meaning but enabling to specify 
a proper form of the potential ip. Here, we demand that for y = expression (fTTj) takes 
the form of Paczyhski-Wiita potential 0, which gives C = 1/2. The PN gravitational 
potential for the SdS spacetime then has the particular simple forrrjj] 



* = - _:' '\ - (12) 



1 + yr 3 /2 
r — 2 — yr 3 

The PN gravitational potential ip diverges at horizons of the SdS spacetime, and 
reaches its maximal value at the static radius of the spacetime where the gradient 

changes its sign from positive to negative values, reflecting both the positions of 
black-hole and cosmological horizons and the changing character of the gravitational 
field which becomes repulsive for r > r s (figure [5]). 

| Another possibility is to demand that at the static radius, where the potential JTTJ has a local 
maximum, ip(r s ) = in analogy with the asymptotic behaviour of the Paczynski-Wiita potential in 
the asymptotically flat Schwarzschild spacetime. This gives C = 1/2(1 - 3J/ 1 / 3 ) implying the PN 
potential in the form 

3 rJ/ l/3 _|_ 2 



2(3j/ 1 / 3 - l)(r - 2 - r 3 y)' 



Such a choice is used in |70| . Both particular formulae of the potential, however, describe the same 
physics. 
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As standard practise in mechanics, we can construct the PN effective potential 
for the radial motion by the relation 

I 2 

v e s(r,0;y,l)=i>(r;y) + — (13) 

2r z sm 9 

where I is the PN angular momentum per particle mass. In the PN approach, as well 
as in the Newtonian theory, the motion of a test particle is confined to a central plane, 
which is conveniently chosen to be the equatorial plane (9 = 7r/2). Circular orbits 
(in the equatorial plane) are again given by local extrema of the effective potential 
v c g (|13j) with = 7r/2. Related energy and angular momentum per particle mass are 
given by the relations 



1 

e c = - 
2 



1 - 



r(r — 3) 



(r — 2 — yr 3 ) 
■^l-yr 3 )} 1 / 2 
r — 2 — yr 



(14) 

, ■ (15) 



2.3. Comparison of general relativistic and pseudo-Newtonian effective potentials 

Despite of the fact that the PN effective potential v c s differs from its GR counterpart 
V c ff , it enables to determine basic features of the test particle motion in accord with 
GR results, especially in the case of the circular orbits. Since behaviour of the angular 
momenta L c and Z c , governing local extrema of the effective potentials V g and v e g, is 
similar, and positions of the L c and l c local extrema coincide, the PN calculation of 
the circular geodesies loci corresponds to the GR results. Marginally stable circular 
geodesies are determined by the condition 

v = MO - r3 (17-15) ( 16 ) 

that is identical with the GR relation [30j . Maximum of the function y ms {r) takes 
the value j/ mSi0 = 12/15 4 , which represents the limiting value for the SdS spacetimes 
admitting stable circular geodesies and, therefore, for the existence of accretion discs 
as well [32] (figure CQ). 

The PN energy per unit mass e c and its GR counterpart E c are related by the 
formula 

-sO-w)- (17) 

Therefore, loci of the marginally bound circular geodesies, corresponding to two 
unstable circular geodesies with the same energy -E c .mb, coincide with two unstable 
circular orbits with the same energy e Ci , n b and angular momentum l c ^ m h — i c .mb/-Kc,mb 
obtained by the PN approach. 

A detailed comparison of the GR and PN approach related to the test particle 
circular motion is presented in [70j and will not be repeated here. We only summarize 
that for y < 10~ 6 the PN formalism reflects quite well the properties of test particle 
motion down to the marginally bound circular orbit with precission growing with 
the cosmological parameter y descending. On the other hand, the PN gravitational 
potential does not reflect the existence of the photon circular geodesic and the 
properties of unstable circular geodesies under the inner marginally bound circular 
geodesic. Nevertheless, these orbits are irrelevant for accretion discs and their 
properties discussed in the present paper. 
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3. Equilibrium configurations of barotropic perfect fluid tori 

In studies of thick accretion discs with non-negligible pressure gradients, investigation 
of equilibrium configurations of perfect fluid tori plays a significant role [72]. The 
tori are considered to be configurations with negligible self-gravity and no influence 
on the spacetime structure, determined by the BH spacetime geometry and an 
appropriately chosen rotational law. We assume a barotropic fluid, i.e. a fluid with the 
pressure-energy density relation of the form p = p(e), for which, according to Boyer 
[75], the stationary equilibrium configurations are described by closed equipotential 
surfaces of the 'gravito-centrifugal' potential. Thick discs are characterized by toroidal 
equipotential surfaces. The necessary condition for the accretion process to be allowed 
is the existence of critical, marginally closed equipotential surface with an inner cusp, 
which enables an outflow of matter from the disc into the black hole. In spacetimes 
with a repulsive cosmological constant (A > 0), the critical point (cusp) of the 
marginally closed equipotential surface can be located also at the outer side, enabling 
thus the outflow of matter from the torus into the outer space. Such a configuration is 
called excretion disc. In the special case, both cusps exist at one critical equipotential 
surface, enabling simultaneous outflows of matter from the disc into the black hole 
and to the outer space [32] [33] . 

In the following, we show that the PN gravitational potential (fl"2|) enables to 
construct barotropic perfect fluid tori in the framework of the PN theorjfjj with 
properties very close to relativistic thick discs. 



3.1. General relativistic perfect fluid tori 

The stress-energy tensor of a perfect fluid with an energy density e and a pressure 
p (as measured by comoving observers) orbiting a black hole is given by the relation 
T lk = (e + p)U" l U k + pg lk . The fluid is assumed to move in azimuthal direction only, 
thus the 4- velocity field of the fluid has only two non-zero components: C/'(r, 6*) and 
U^(r,6). Projecting the covariant energy-momentum conservation law VfeT lfc = 
onto the hypersurface orthogonal to the 4- velocity U\ we obtain the relativistic Euler 
equation in the form [75j [76] 



=-diQnU t ) + ^L, l = {r,8} (18) 



p + e n tJ l-Slf 
where in general stationary and axisymmetric spacetime 

t z gu + 2lg t </, + g<f,<p 

The angular velocity (related to distant observers) = U^/U 1 and the specific angular 
momentum i = L/E = —U^jUt are related by the formula 

n = _ £gu + 9t^ (2Q) 

f-SH + 94><f> 

In the SdS spacetimes, relations (fTT)]) and (|20p reduce to the form 

(tt \2 _ 9tt94>4> ^ 9tt , 01 n 
\ u t) — -75 ; , -j — ■ l/JJ 

« 9tt + 94»i> * 9<P<P 



% Within the PN approach, the thick accretion discs orbiting the Schwarzschild black hole were 
treated in detail, e.g., in | 74| . 
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For any barotropic fluid, i.e., the fluid for which p = p(e), a solution of the 
relativistic Euler equation can be given in terms of Boyer's condition for constant- 
pressure surfaces, which coincide with equipotential surfaces of the potential W(r, 9) 
defined by the relations [75] ES] 

r p dp 



P 



W in - W, (22) 



W- m -W = \n (U t ) in - In (U t ) + jf ^-^L; (23) 

the subscript 'in' refers to the inner edge of the disc in the equatorial plane. The 
equipotential surfaces are determined by the condition W{r, 9) = const and in a given 
spacetime can be found from the equation (|23p . if a rotational law O = is given. 

In the case of discs with uniform distribution of the specific angular momentum 
£(r, 9) = const0, the potential W is given by a simple formula 

W(r, 9) = In U t = ln t J r - 2 - yr T r ZL, - (24) 



\r 3 sin 2 



( r _ 2-yr 3 )e 2 ] 1 / 2 



The first reality condition, r — 2 — yr 3 > 0, restricts the existence of equipotential 
surfaces to stationary regions of the spacetime. The second reality condition, 
r 3 sin 2 9 — (r — 2 — yr 3 )£ 2 > 0, evaluated in the equatorial plane (9 — n/2), implies a 
restriction given by the photon motion 

£ 2 <e 2 ph (r ]y )^ f 3, (25) 

as the function ^ 2 h (r; y) plays the role of the effective potential governing the photon 
geodesic motion in the equatorial plane [3D]. The function £p h (r;y) has one local 
extreme (minimum), ^ 2 hc , located at r = 3 independently of the cosmological 
parameter y, and corresponding to the circular photon orbit in the equatorial plane. 

Character of the equipotential surfaces is fully determined by the existence of 
local extrema of the function W(r, 9) and their mutual behaviour. Since the orbits 
with vanishing gradient of W, i.e., satisfying conditions d r W(r,6) = d$W(r,6) = 0, 
correspond to places with zero pressure gradients, the fluid has to follow the geodesic 
motion there. As shown in [32], the local extrema of W are located only in the 
equatorial plane (9 — tt/2), where the radial profile of the potential is described by 
the function 

W 7r/2 (r) = W(r,e = n/2). (26) 

The necessary condition for the local extrema, d r W 7r /2( r ) = 0, implies 

_ r 3 (l — yr 3 ) 
(r — 2 — yr 3 ) 2 

Since £ 2 (r;y) corresponds to the specific angular momentum distribution of circular 
geodesies, the extrema of the potential (corresponding to circular geodesies) fulfill the 
relation 

W ext = W„ /2 (r; £=£ c ,y) = \nE c (r; y), (28) 
where E c is the energy per unit rest mass related to the circular geodesies (J7J. 

+ This case corresponds to the marginally stable tori with properties reflecting quite well those of 
tori with a general allowed specific angular momentum profile. 



e=ll{r;y)= r l ^ " J . (27) 
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Figure 3. Functions (r; y) (solid) and £^ h (r; y) (dashed). There are four classes 
characterizing mutual behaviour of these functions (Al, A2, A3, B) and two 
different types of behaviour of the function £^(r;y) (a,b). The horizontal dashed 
and dotted lines denote values of and £^ h c for a given cosmological parameter 
y, respectively. 



There are four types of the £ 2 h and l\ mutual behaviour, depending on the value 
of the cosmological parameter y (figure [3] and table [T]) . They give four classes of 
the SdS backgrounds, differing in the types of the structure of equipotential surfaces 
admitted by the spacetime. For y ms . c < y < y cr n = 1/27 = 0.037, the function 
£^(r; y) is monotonically decreasing, corresponding to unstable circular geodesies only. 
For y < y mSje = 12/15 4 = 2.4 x 10 -4 , the function £ 2 (r;y) has two local extrema, 
the minimum ; and the maximum £ 2 as a , corresponding to the inner and outer 
marginally stable circular geodesies. The rising part of £ 2 (r;y), located between the 
local extrema, corresponds to stable circular geodesies, while the descending parts of 
£ 2 c {r;y) to unstable ones. If y = y e = 1/118125 = 8.5 x 10~ 6 , £ 2 ms o = ^ h c , and for 

V < Vei £ ms .o > ^ph,c- 

Equilibrium toroidal configurations of a barotropic fluid can exist only in the 
SdS spacetimes enabling motion along stable circular geodesies, i.e. the spacetimes 
with y < y m s,e, and with the specific angular momentum restricted to the range 
£ G (4ns, i, ^ms,o)- For such a constant specific angular momentum distribution, the 
radii given by the condition £ = £ c (r; y) correspond to the motion of the fluid along 
unstable, stable and again unstable circular geodesic, as the radius grows. The 
stable circular geodesic represents a center of the torus, because the potential W(r, 9) 
has a local minimum there, while the unstable circular geodesies determine critical 
points (cusps) where the potential W{r, 9) has local maxima and the corresponding 
equipotential surface is self-crossing. If the critical surface with the cusp is marginally 
closed, i.e. it encloses other closed equipotential surfaces corresponding to lower 
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Table 1. Classification of the SdS spacetimes. There are four types of the 
SdS spacetimes differing in the behaviour of the functions t^.(r\y) and £^ h (r;y) 
(relevant for the GR classification) and two types of behaviour of function l^{r; y) 
(relevant for the PN classification). Note that in realistic astrophysical situations 
the case Al is relevant. 



PN 


Relativistic 


Range of y 


class 


class 




a 


Al 


0<y<y c 


a 


A2 


V = Ve 


a 


A3 


Ve < V < S/ms.c 


b 


B 


Vms.c <V< 2/crit 



values of the potential W(r,6), an outflow of matter from the torus through the 
cusp is possible. For £ € (Cs.i,^mb) where £ m b is the specific angular momentum of a 
particle moving along the inner/outer marginally bound circular geodesic, the critical 
marginally closed equipotential surface self-crosses in the inner cusp; this configuration 
corresponds to the accretion disc. For I e (Cb^ms,o), the critical marginally closed 
equipotential surface self-crosses in the outer cusp; this configuration corresponds to 
the excretion disc. For £ = £ m b, the critical marginally closed equipotential surface 
possesses both cusps, therefore, common accretion onto the BH and excretion into 
the outer space are equally possible; this configuration is called 'marginally bound 
accretion disc'. When y < y c and £ p h :C < £ < £ mS:0 , a forbidden region for the 
fluid with a given £ — const exists, restricted by the radii following from the relation 
£ = £ pb (r;y). 

Due to the properties of the functions £^, h {r;y) and £l(r;y), we can deduce 
that there exist twelve characteristic structures of the equipotential surfaces; for the 
astrophysically most relevant relativistic class Al (see table [I]), nine characteristic 
structures of the equipotential surfaces exist and five of them correspond to toroidal 
configurations (for details see [32] )• In dependence on £, the configuration describes the 
accretion disc, marginally bound accretion disc, excretion disc without the forbidden 
region, excretion disc with the forbidden region in the form of infinitcsimally thin 
shell, and the excretion disc with the finite forbidden region. 



3. 2. Pseudo- Newtonian perfect fluid tori 

The PN version of the Euler equation for the perfect fluid characterized by the pressure 
p and the mass-density p, rotating in a central gravitational field, takes the form [69] 

-d iP =-d iW , i = {r,6}, (29) 
P 

where the potential w is defined by the relation 

r i 2 

w = ^j-Jj^dR, (30) 

and R = rsin#. Again, for a barotropic fluid, p = p(p), the surfaces of constant 
pressure are just the equipotential surfaces of the potential w: 

rPd P *m 
— =w in - w. (31) 

P 
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In the SdS spacetimes, assuming the PN gravitational potential ip given by relation 
(|12p and uniform distribution of the angular momentum per unit mass, I = const, we 
obtain the potential w in the for n0 

K ' r-2-yr 3 2r 2 sin 2 9 ' 



which is equal to the PN effective potential governing the radial motion (JTJ 

The PN potential w diverges at horizons of the SdS spacetime and is regular in 
the whole stationary region of the spacetime, but it does not reflect the existence of 
the photon circular orbit [70) . Its radial profile in the equatorial plane is described by 
the function 

w v/2 (r) = w(r,9 = w/2). (33) 

Extrema of w„/2 (r) correspond to the motion along Keplerian circular orbits with the 
energy per unit mass e = e c and the angular momentum per unit mass I — l c given 
by HID and (S), 

w cxt = w 7T / 2 (r; I = l c , y) = e c (r; y). (34) 

In the case of the PN potential w, the structure of its equipotential surfaces 
is governed by the function i\(r\y) only, and we can find seven various structures 
of equipotential surfaces in dependence on the parameters y and I. For y < y ms , e 
and I € (ims,i,ims,a), the closed equipotential surfaces of the potential w exist, 
corresponding to barotropic fluid tori. The center and cusps are given in exactly the 
same way as in the GR framework. Again, we obtain configurations corresponding to 
accretion discs, marginally bound accretion discs and excretion discs (figure 2]). 

3. 3. Comparison 

Topology of equipotential surfaces of the potentials W and w can be directly inferred 
from their behaviour in the equatorial plane, W w /2(t] y) and w„/2{r; y). Note that the 
local extrema of both the functions are governed by the same function £ c (r; y). Perfect 
fluid follows circular geodesies along the radii where the conditions I = l c (r;y) and 
£ — £c(r; y) are satisfied, which is obvious from the fact that the conditions d r W n /2 = 
and d r w n /2 — are equivalent to the conditions of vanishing of pressure gradients. 

Our construction of the PN gravitational potential ip guarantees that the function 
£ 2 (r; y) describes the specific angular momentum of freely moving particles in both the 
GR and PN approach. Therefore, we can expect close similarities of the character of 
equipotential surfaces constructed in both approaches. It can be even shown that the 
shapes of the equipotential surfaces constructed in GR and PN approaches coincide 
(figure 2]). Really, the equipotential surfaces W(r,6) = const and w(r,9) = const are 
determined by the same relation 

d0 = _drW _8rW 

dr 8gW 8gW [ ' 

_ r 3 (l - yr 3 ) sin 2 9 - (r - 2 - yr^ 2 " 2 
r(r — 2 — yr 3 ) 2 £ 2 

* Constant of integration in J30D is chosen to be zero. 



tan( 



(loer) s'm8 




-2-10 1 2 



(logy) sin# 

(b) 

Figure 4. Meridional sections of equipotcntial surfaces (potential contours) of 
the potentials W(r, 6) and w(r,9), and their behaviour in the equatorial plane 
(8 = it/2) represented by the functions W^^M^ashed curve and /2( r )~ sou d 
curve. Self-crossing (critical) equipotcntial surfaces correspond to the values of 
potentials at their local maxima, (a) £ ms ; < £ < £ rn \ 1 : the inner critical surface is 
marginally closed, the outer critical surface is open, (b) £ = £ m b : the only critical 
surface containing both cusps is marginally closed, (c) £ m t, < I < <? p h, c < "ms,o 
or Anb < (■ < 'ms,o < ^ P h,c : the inner critical surface is open, the outer critical 
surface is marginally closed, (d) £ p ^ tC < I < < mB)0 : In GR approach, there 
is a 'forbidden region' for the fluid with prescribed specific angular momentum 
and, thus, no inner (open) critical surface; the outer critical surface is marginally 
closed. In PN approach, the 'forbidden region' does not exist, the inner critical 
surface is open, the outer critical surface is marginally closed. 
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Figure 5. Potential barriers of (a) relativistic thick discs in dependence on the 
specific angular momentum £ £ (^ms(i) > ^ms(o) ) an d (b) pseudo-Newtonian thick 
discs in dependence on the angular momentum per unit mass I £ (i ms (i) , ' ma (o) ) f° r 
three fixed values of the cosmological parameter y: 10 — 6 (dotted), 10~ 5 (dashed) 
and 10~ 4 (solid). Accretion/excretion discs are represented by thin/thick curves, 
the maxima correspond to marginally bound accretion discs with 1 = 1 = £ m h . 



following from the conditions dW = and dw = 0. Using the Schwarzschild 
coordinates in GR approach and spherical coordinates in PN approach, the 
equipotential surfaces have the same shape determined by the relation 

-l 



sin 2 I 



(36) 



C ~ z?2 (21 Tf. _o_„,z?3V ( 37 ) 



l 2 (r- 2-yr 3 ) 
where the integration constant can be expressed as 

1 Ru 

i?.f n e(B, in -2-yRl 

In the PN approach we assume a flat space, while in the GR approach, we have to 
take the spacetime curvature into acount. Therefore, the proper space distances inside 
the PN and GR tori differ and we can expect distinct potential barriers (potential 
differences in the center and in the edge of the torus) or different total masses of the 
torus, as predicted by the GR and PN approaches for otherwise identical conditions. 

The potential depths comparison can be realized in the following way: 
determining values of the potentials W and w at their extrema in the equatorial 
plane, i.e., their inner maxima, W e xt,i and w ex t,i, located at the inner cusp of 
equipotential surfaces, minima, W cx t,c and u> ox t.c, located at the central ring of 
equipotential surfaces, and outer maxima, W ex t,o and w cx t,o, located at the outer 
cusp of equipotential surfaces, we can determine the related GR and PN inner and 
outer potential barriers (figure [5]) 

Ui - [W ext ,i(^ y) - W eKtiC (l; y)} , (38) 

Mi = [w ext< i(l; y) - i/w,c(2; y)] , (39) 

which represent accretion discs, 

U = [W e ^, (£; y) - W c ^ c (£; y)} , (40) 
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Figure 6. (a) absolute and (b) relative differences of GR and PN potential 
barriers in accretion (thin) and excretion (thick) discs for three fixed values of the 
cosmological parameter y: 10 -6 (dotted), 10 -5 (dashed) and 10 — 4 (solid). 



u = [w e xt,o(ky) - Wext, c (hy)] , (4i) 

which represent excretion discs, and their differences 

A; = Ui(£;y) - Ui(l;y), (42) 

A = U (£;y)-u (l;y), (43) 

characterizing variances in the GR and PN descriptions of the accretion and excretion 
discs. Although the maximal absolute differences |A; >0 | grow with the cosmological 
parameter descending, the relative (percentual) differences Ai, Q = |Ai )0 /E/i j0 | x 100% 
behave in opposite way (figure [6]) . 



4. Adiabatic perfect fluid tori 

Important physical information is reflected by the behaviour of the fluid inside the 
toroidal configurations. In realistic accretion discs, it is a very complex problem how 
to treat all relevant phenomena concerning properties of the fluid and generation of 
heat and radiation. Here, in order to compare the GR and PN approaches, we restrict 
our attention to a very simple case of adiabatic discs (with negligible influence of 
radiation) . 

We assume a perfect fluid obeying an adiabatic pressure-mass-density relation 

p = Kp 1+1/n , (44) 

where K and n are the adiabatic constant and adiabatic index, respectively. Further, 
we assume that the pressure is given by the relation for an ideal gas 

P = —k B T (45) 
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where T is the thermodynamic temperature, and constants p., m u and fee are mean 
molecular weight, atomic mass unit and the Boltzmann constant, respectively. In this 
case (no radiation), the energy density is given by the relation 



e = p + np. 
4--1- General relativistic tori 



(46) 



Solving equation (f2"2")l for the adiabatic fluid, we can express distribution of all 
basic thermodynamic quantities throughout the toroidal configuration in terms of 
the potential W in the following way: 
r ^aw _ 

AW = Wi n - W, (47) 



K{n+l) 



K~ 



.A1V 



-1 



Tt + 1 



e =p 



T = 



1 



e" 1 ' - 1) 



r e AW -i - 


n - 


K(n + 1)_ 





,AW 



fim u 



Kp 



l/n _ 



/im u e 



AW 



1 



fee n + 1 

In adiabatic tori, the equipotential surfaces W(r, 9) 



(48) 
(49) 
(50) 
(51) 

const also determine (and are 



identical) the surfaces of constant thermodynamic quantities p, e, p and T. Extremal 
values of the thermodynamic quantities are reached at the centre of the tori. 

For a given mean molecular weight /i, which depends on the matter of the disc, 
relation (|51[) enables to express the temperature in the centre of the torus in the form 



T c , = 



p,m u e 



-1 



AW r = W m -W c . 



fee 7i + l 

Temperature profile of the torus is thus given by 
T = T, 



e AW c _i 

Specifying the central mass-density p ci the adiabatic constant K, the central 
energy density e c and the central pressure p c are given by the relations 



(52) 



(53) 



A 

Pc 



-1/r 



1 



n + 1 

71 



n + 
e A ^-l 



>AW C 



(54) 
(55) 

Pc=Pc — — • (56) 

n + 1 

In terms of p c , which is the only free parameter in description of adiabatic tori, the 
mass-density profile is given by 

/ e AW _j \ » 

P = Pc[ zxw—r ) ■ ( 57 ) 



gAWc _l 
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Similarly, the pressure profile is described by the relation 



P=Pc 



5 AW 



-1 



e AW c _! 



(58) 



Central values of the thermodynamic quantities are given by the potential 
difference between the edge and the centre of the torus, A c = W- m — W c , which grows 
with the £ = const growing from £ m s,i to £ ms ,o- For £ = £ m b, the toroidal structure 
is the most extended one and the central values of the thermodynamic quantities 
reach their extreme values in a fixed BH background. Analogical statements as in GR 
approach hold, of course, in the PN approach. 



4-2. Pseudo-Newtonian tori 



v -- 

T 



K- 



K(n + 1) 

Win 



71 + 1 

fim u w m — w 



n+1 



Starting with relation (|3ip . we can express distribution of all basic thermodynamic 
quantities throughout the adiabatic torus in terms of the PN potential w in the 
following way: 

-i n 

" (59) 

(60) 

k B n+1 ■ < 61 > 

Again, the equipotential surfaces w(r, 9) — const are identical with surfaces of constant 
thermodynamic quantities that reach maximal values at the centre of the torus. 
The temperature profile of the torus is described by the relation 

Win - W 



T = T C 



(62) 



Win - W c 

where the central temperature of the torus is given by 

Hm u Win - w c 

Kb n+1 

In terms of the mass-density in the centre p c , the adiabatic constant K, the central 
pressure p c and the mass-density profile are given by the relations 



K = p- 1/n - 



n -\ 

- Wt 



1 



w 



Pc , 

s Win - W, 

The pressure profile is then given by 

Win - w 



P=Pc 



n+1 



(64) 
(65) 

(66) 
(67) 
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Figure 7. PN (solid) and GR (dashed) temperature profiles of adiabatic 
marginally bound accretion discs in the equatorial plane. The tables give a GR 
and PN temperature in the centre of the disc for a given adiabatic index n. 

4-3. Comparison 

In order to compare the PN and GR descriptions of adiabatic tori, we plot the 
temperature and mass-density profiles of the GR and PN marginally bound accretion 
discs (tori with £(r,9) = £ m h) in the equatorial plane (6 = tt/2) for various values of 
the adiabatic index n and cosmological parameter y (figures [7] and [8]). Maximal values 
of T and p correspond to the centre of the torus. We can see that the PN-temperature 
profile is little bit higher than the GR-temperature profile and the difference of these 
values at a given r grows with growing cosmological parameter y. On the other 
hand, the PN density profile p/p c is only slightly lower than the GR density profile, 
and the difference of these values descends with the cosmological parameter growing. 
The difference is demonstrated by the quantity Aqr-pn in figure [3 It should be 
noted, however, that when the adiabatic constant K is the same for both PN and GR 
discs, the mass-density predicted by the PN approach (relation ([59]) ) is always higher 
than the mass-density given by the GR approach (relation (|47|) ) at a given r. Their 
difference is the biggest in the centre of the torus. Ratio of the PN and GR predictions 
of mass densities in adiabatic tori is shown in figure [9l 

Similarly, we can also plot pressure profiles of GR and PN adiabatic discs for 
various values of the adiabatic index n and cosmological parameter y, described by 
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Figure 8. Equatorial PN and GR mass-density profiles of adiabatic marginally 
bound accretion discs in units of the central density p c . Difference of the GR and 
PN values of p/p c (denoted as Aqr_pn) is shown by small figures for various 
values of the adiabatic index n. 



relations (|58[) and (|67p . For given values of y and n, the profiles p/p c show the 
same behaviour as mass-density profiles for the same cosmological parameter y but 
adiabatic index (n — 1). Thus, the pressure profiles are narrower than corresponding 
mass-density profiles. 

In the framework of general relativity, it is possible to analyze how much 
relativistic is the fluid matter in the adiabatic torus in dependence on the adiabatic 
index n and cosmological parameter y, investigating relation (|49p for the special energy 
density. The results are shown in figure [TOl We can see that the importance of the 
GR effects in the torus grows with the adiabatic index n growing and the cosmological 
parameter y descending. However, the special- relativity effect is important only in 
the inner part of the disc around its centre, and the total energy density in the centre, 
where it is maximal, exceeds the rest-energy density maximally by ~ 3% of the rest 
value. 

5. Total mass of adiabatic tori 

In this section, we compare GR and PN results giving a total mass of the adiabatic 
perfect fluid torus. Using this quantity, we can put some limit on validity of the 
test-disc approximation implicitly assumed in our calculations. 
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5.1. General relativistic mass formula 

Total mass of a torus is determined by the Tolman formula [75] 

m = [ (-2? + I7 + I? + 2j)V=?drd0d0 ) (68) 

J disc 

where are components of the perfect fluid stress-energy tensor and g = det(gik). 
For a perfect fluid torus orbiting a SdS black hole, the Tolman formula reads 



disc 



l-ne 



(e+p) + 2p 



• sin 9 drd9dcj>, 



(69) 



where the angular velocity profile f2(r, 0) is determined by the spacctimc geometry 
and the specific angular momentum distribution £(r, 9) through the relation (I2ip : 

r-3 



fi(r, 9) 



r-2 



r 3 sin 2 i 



l{r, 9). 



(70) 
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In the case of the adiabatic tori characterized by the adiabatic index n and central 
density p c , relation (|69p can be rewritten to the form 



Pc 



disc . 

AW i 



1 - Hi n + 1 



X {^Aw7^l) r 2 smedrd9d^ (71) 

enabling to find the mass of the adiabatic torus from the behaviour of the potential 
field W(r, 9) in the disc. In the limit of non-relativistic fluid for which p <C e w p, the 
mass of the adiabatic torus is approximatively given by 

_ r i + nt ( w in - w \ n a . 

m ~Pc q ^ I 7J7 fjT" ) r sin 6> drd6>d</>. (72) 

•/ disc 



1 - qm w in - w c 



5.U. Pseudo-Newtonian mass formula 

In the Newtonian description of the fluid, we naturally assume a non-relativistic perfect 
fluid with p <C p- Geometry of the spacetime is Euclidean. 
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The mass of a torus can be determined from the relation 
m = p(r, 9, 4>)r 2 sin 9 drd9d(f>. 



(73) 



disc 



Using the results of the previous section for adiabatic PN tori, the mass of such a torus 
is given by the formula 



5.3. Comparison 

Information on mass of the torus provides an opportunity to discuss the relevance of 
the test-disc approximation. It is not trivial to find a reliable criterion in this case 
because of large extension of the torus in comparison with the BH extension. A very 
simple criterion is to compare the total mass of the torus with mass parameter of 
the black hole, stating that the test-disc approximation is acceptable for m <C M . 
Due to much smaller 'compactness' of the torus, as compared with the black hole, 
we assume m < M/10 to be sufficient enough for the test-disc approximation, taking 
m = M/10 as a testing limit. Central mass densities of adiabatic tori with m = M/10 
are presented in table [2j Thus the test-disc approximation should be relevant for 
adiabatic tori with central mass densities lower than those in the tables (for given 
adiabatic index n and mass of the black hole M). 

In order to compare the GR and PN masses of adiabatic tori, we calculate in 
both the approaches the central mass densities, corresponding to the same mass of 
the torus. At first, we choose the marginally bound accretion discs (which correspond 
to the most extended toroidal structures in any given spacetime), i.e., the discs with 
L = ^mb, and calculate the mass densities in the centre for the limiting (test-disc) 
situation when the total mass of the torus is given by 10m = M = y3?//A = 1 in 
dimcnsionlcss formulation. 

With respect to the current value of the cosmological constant (Ao = 1.3 x 
10 -56 cm -2 ) and recent observations confirming existence of black holes with masses 
up to M ~ 1O 1O M0 [66], the highest current astrophysically relevant values of the 
cosmological parameter are about y ~ 10 -25 . For marginally bound accretion tori 
orbiting intermediate-mass and supermassive black holes, the above-mentioned central 
mass-densities (and pressures) are presented in table [5] in physical units. Ratio of such 
PN and GR central mass-densities in dependence on the cosmological parameter y is 
presented in figure [Til for three typical values of the adiabatic index n. Notice that for 
astrophysically relevant black holes with y < 10~ 25 the precision of the PN approach 
is relatively high and differences of GR and PN results are few percent only. 

Analogical calculations were performed for accretion discs with £ = l/2(£ ms ; + 
£ m b) and excretion discs with I = l/2(4nb + ^ms,o)- The results are again presented 
in table H 

We can see that the PN densities and pressures predominate over the GR ones, 
but both the approaches imply results of the same order. In the case of excretion 
discs, the results even almost coincide. Apparently, this is because the excretion discs 
do not take place in 'highly-curved' region of the spacetime nearby the black hole. 
For maximally extended tori (£ = £ m b), differences of the PN and GR approaches 
are relatively small (few percent), since only small part of such tori are located in the 
'highly-curved' region close to the BH horizon. The PN and GR results are of the same 




(74) 



£ 

Table 2. GR and PN central mass densities and pressures of adiabatic thick discs with total mass m = M/10; M is the mass parameter of o 

a particular SdS black hole. The results are given for current value of the cosmological constant, Ao = 1.3 X 10 — 56 cm" 2 , and three values |> 

of the adiabatic index n = (1.5, 2.5, 3.0). Size of the disc is characterized by position of the outer cusp of toroidal configuration i? ut (given g 

in units of M and in parsecs). The inner cusp of accreting configurations is located in vicinity of the black hole at R ln 4M, while the g" 

inner edge of excretion discs is located far away from the black hole at ~ 10 7 M for y = 10~ 26 , or ~ 10 10 M for y = 10~ 34 . § 
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Figure 11. Ratio of central mass densities of the PN and GR adiabatic marginally 
bound accretion discs of the same mass equal to the mass of the black hole, 
i.e. m = M, for three values of the adiabatic index n = (1.5, 2.5, 3.0) as a 
function of the cosmological parameter y. Attached table shows the ratios for 
three representative values of y. 



order (tens of percent differences) even for accretion tori fully located in the strong 
gravity region nearby the horizon. With the cosmological parameter y descending, the 
precision of PN calculations increases. 

Comparing results in table [2l obtained for the case m — M/10, with relations 
(jTTj) and l[71j) for mass of the torus in both approaches, we can deduce that in 
the case of equal central mass densities of PN and GR tori, the mass of PN torus 
must be a little bit smaller than the mass of its GR counterpart (but of the same 
order!)0 For example, if we choose marginally bound adiabatic torus with n = 1.5 
and p c = 10~ 15 kg/m 3 , orbiting a supermassive black hole with M = 10 6 M©, the 
mass of such a torus in the PN approach will be topn = 2.8 x 10 -5 M, while in the 
GR approach ttt-gr = 2.9 x 1CP 5 M. Again, with decreasing values of the cosmological 
parameter y, precision of the PN calculation increases. 

6. Discussion 

We shall discuss behaviour of thermodynamic quantities in three types of toroidal 
configurations, and present some comments on the definition of the PN potential. 

6.1. Profiles of thermodynamic quantities 

Central values of mass-density p, pressure p and temperature T in adiabatic tori, 
suggested by PN and GR approaches (tables [HO, are given for current astrophysically 
relevant values of cosmological parameter y E (1CP 44 — 10 -26 ), reflecting the range 
of BH masses from stellar-mass (10 M Q ) to the most supermassive ones (1O 1O M , 
which is comparable with mass of the heaviest supermassive black hole TON 618 

tt Since the PN approach is denned in order to describe the BH physics in the Euclidean space, the 
proper volume of the curved region in the GR approach is greater than its flat PN counterpart. 
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Table 3. GR and PN central temperatures of adiabatic thick discs with total mass 
m = Af/10; M is the mass parameter of a particular SdS black hole. The results 
are given for current value of the cosmological constant, Ao = 1.3 X 10 — 56 cm -2 , 
and three values of the adiabatic index n = (1.5, 2.5, 3.0). 
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being observed till now) and the present value of the cosmological constant Ao = 
1.3 x 10~ 56 cm~ 2 . The calculations are performed for both accretion and excretion 
discs, as well as for the most extended marginally bound tori enabling simultaneous 
outflow of matter into the black hole and to the outer space. In all cases we assume 
the test-disc approximation, with the limiting value m = M/XO. 

6.1.1. Most extended tori In the marginally bound tori with i = const = £ m b 
orbiting any astrophysically relevant SdS black hole, the central temperature reaches 
the extreme value independently of the mass of the torus, since for adiabatic fluid it is 
fully determined by the potential depths AW C = Wi n — W c in GR and Aw c = w m — w c 
in PN approaches that reach their maxima for the tori with I = £ m b (cf. equations 
(|52[) and (frJ3")) ). Note that for the mentioned interval of y the central temperature is 
almost independent of y for both GR and PN models, being a few percent higher in 
the PN model as compared to GR ones. On the other hand, the central values of 
mass-density and pressure depend on the total mass of the disc (cf. equations ([7l"|) . 
(j74|) and ([56]). ([65]) ). Extension of I = £ m h tori expressed in terms of gravitational 
radius of the central black hole M falls with the cosmological parameter y growing. 
Since for y < 1CP 10 , there is r m b i0 ~ r B = y _1 ^ 3 M and r^ ~ M, an approximate 
relation for r m b, can be given in a simple form r m b i0 /ni ~ M -2 / 3 A~ 1 / 3 . The central 
mass-density and pressure grow sharply with descending y (BH mass). The PN model 
predicts the central density and pressure by < 10% (and typically only by few percent) 
higher than the GR model. 

6.1.2. Accretion tori For £ ms ,i < £ — const < £ m b, both the extension of the accretion 
disc (toroidal structure with the inner cusp) and the central potential depths AW C (or 
Aui c ) are smaller as compared with the most extended marginally bound accretion 
discs. Thus the central temperature is also smaller. On the other hand, when the 
condition m = M/10 is assumed, the mass-density and pressure of the accretion tori 
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Table 4. Position of the outer edge i? ut of accretion tori in units of M in 
dependence on the value of constant specific angular momentum (the results are 
the same for both PN and GR marginally stable tori) and the GR values (of the 
same order as PN values) of central densities in kg/m 3 for n = 1.5 in the disc which 
is characterized by parameter k in the relation t = ^ mSl i + fc(<? m b -f m s,i) and value 
of cosmological parameter y = 10 — 26 . Note that there is no significant dependence 
of the inner edge position on the parameter k, being located at i?; n = 4M. 
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become higher as compared to the most extended discs. The extension of the toroidal 
structure falls sharply with a small descend of the £ — const value from £ mrj . On the 
other hand, it falls slowly while £ = const approaches £ ms .i, as shown in the case of 
the most massive black hole in the table [4] This is caused by the behaviour of the 
potential W x / 2 (r) which contains an extremely extended and flat plato near the static 
radius r s and deep valley near the inner cusp f figure [T2"]) . When the specific angular 
momentum £ = const — > £ mSi i, the extension of accretion torus shrinks to zero being 
concentrated at r ms .j. 

It is useful to show behaviour of central values of the mass-density and pressure 
in the case of accretion tori restricted to vicinity of the central black hole, i.e., with 
£ = const close to ^ m s,i- In fact, such tori are obtained with a moderate case of the 
representative value of £ £ (£ ms ,i, £mb), namely £ = l/2(^ mS a + £ m b), see table[2] Here, 
we have included also the astrophysically very important case of stellar mass black 
holes, having a typical mass M — 10 M Q , as frequently observed in Galactic X-ray 
binaries [77] . 

Notice that for such values of y and £, the extension of the tori in units of BH 
mass M is almost independent of the cosmological parameter y. However, the central 
mass-density and pressure depend strongly on y and only slightly on the adiabatic 
index n. Contrary to the case of most extended tori, the accretion tori located close 
to the BH horizon reflect strong relativistic phenomena. Therefore the predictions of 
the PN and GR models differ for more than 10% and up to 40%. The chosen value 
of £ = l/2(£ mSi i + £ m b) = const corresponds to tori located close to the BH horizon 
in regions where the influence of the cosmological constant is very small; extension 
of these tori is again small (i? ut = 20.6 M). Since we use the limiting mass of the 
disc m = M/10, central mass-density and pressure reach high values, which grow 
enormously with lowering mass of the central black hole. For M = IOMq, i.e., a 
typical stellar mass black hole, the central density of the torus reaches the values for 
which the neutrons start to drip out of the nuclei [75] . 

6.1.3. Excretion tori For £ m b < £ — const < £ m s,o> both the extension of the excretion 
tori (toroidal structures with the outer cusp) and their central potential depths are 
again lowered in comparison to the most extended tori for fixed BH mass M, and the 
central temperature does not reach values of marginally bound discs. Contrary to the 
case of accretion tori, the central values of mass-density and pressure of the excretion 
tori are also lowered relative to the most extended tori. Extension of excretion tori 
falls very sharply when £ = const slightly overcomes £ m b (due to radical shift of the 
inner edge), while it slowly shrinks to zero extension at r ms o when £ = const — » £ mSt0 
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Figure 12. Behaviour of the relativistic potential W in the equatorial plane of the 
SdS spacetime with y = 10 — 26 for two values of the specific angular momentum: 
£ = £ mb = 4M (solid) and I = l/2(^ mSji +l mb ) = 3.8 M (dashed). Notice the 
extremely fiat plato extended around static radius r s ~ j/ -1 / 3 . Extension of such 
a plato grows with y descending. 



(see table [5] constructed for the most massive black hole with y = 10~ 26 ). Again, 
this is caused by the very flat and extended plato of the potential W v / 2 {r) near the 
static radius (see figure fT3|). While i — const grows for excretion tori of fixed mass 
m = M/10, radius of their centre in the symmetry plane grows, while their section 
perpendicular to the symmetry plane shrinks. The central density falls with I = const 
growing for large part of the range of £ S (Cb^ms.o), but it starts to grow again as 
£ — const is close enough to £ ms ,o, reaching high values especially for £ — > £ ms , , since 
the perpendicular section of the tori shrinks to zero in such a situation (see table [5]). 

We demonstrate behaviour of the mass-density and pressure for excretion tori with 
a moderate value of I = const = l/2(£ m b+£ ms , ). Such a value of £ = const corresponds 
to a torus with rather extreme properties; notice the values of the critical angular 
momentum for y = 1(T 26 , i.e., Cs,i = 3.7 M, £ mh = 4 M and £ ma , Q = 1.5 x 10 4 M. 
The results are given in table [2] The excretion tori are discussed in the astrophysically 
relevant cases of supermassive black holes, when extension of such a toroidal structure 
is comparable to the related galaxy extension lffl 

The central mass-density in such excretion tori is in many orders lower than those 
of the most extended tori (or the accretion tori) having the same mass m = M/10. It 
is almost independent of the cosmological parameter y, and thus on the BH mass, 
contrary to the other two cases, where the differences are in many orders. The 
excretion tori with such a moderate value of i = const are located in regions close to the 
static region, where the PN approach can be applied with an extremely high precision 
- in fact, the PN and GR predictions differ at the 7-th (or higher) decimal order. 
Further, the central densities and pressure are almost independent of the adiabatic 
index n for such distant excretion tori. 

ff Recall that for intermediate or stellar-mass black holes in star clusters or binary systems, the 
mentioned toroidal configurations are much more extended in comparison with the clusters or binary 
systems. 
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Table 5. Position of the inner edge R ln of excretion tori in units of M in 
dependence on the value of constant specific angular momentum (the results are 
the same for both PN and GR marginally stable tori) and the GR values (of the 
same order as PN values) of central densities in kg/m 3 for n = 1.5 in the disc which 
is characterized by parameter k in the relation I = <? m s,o — fc(£ m s,o — t-mh) and value 
of cosmological parameter y = 10~ 26 . Note that there is no significant dependence 
of the outer edge position on the parameter k, being located at R ou t ~ 10 s M. 
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6.2. On definition of pseudo-Newtonian potential 

Some astrophysical applications of exact GR methods represent samples of too 
complicated calculations. Therefore, approximative PN or even Newtonian methods 
are usually used in order to obtain more or less satisfactory results. There is a variety of 
PN approaches based on different definitions of the gravitational potential, describing 
some aspects of BH physics [69j EU EH EH ED [82] ■ 

In study of accretion discs, the most convenient and widely used is the Paczyriski- 
Wiita gravitational potential |69j . Originally, this potential was introduced by a guess, 
when attempting to include the Schwarzschild radius into the Newtonian gravity. 
However, there is a simple heuristic method for derivation of the PN potential based on 
the properties of the test-particle motion, giving in the case of Schwarzschild spacetime 
just the Paczyhski-Wiita potential and used in derivation of the PN gravitational 
potential for the equatorial plane of the Kerr spacetime as well [71] . In this paper, we 
used this method to define the PN potential ip for the SdS spacetimes. 

Looking for a definition of PN potential from the point of view of fluid tori 
description in spherically symmetric spacetimes, one can state the following. The 
Newtonian and GR versions of potentials, determining equipressure surfaces in 
barotropic marginally stable tori, are given by 

J! 

2r 2 sin 2 



WN = -0N + — 2 - , (75) 



W = ln( ~ 9a9 ^ Y /2 , (76) 



where V'n and I are the Newtonian gravitational potential and angular momentum 
per unit mass, while I = L/E is the relativistic specific angular momentum. For 
marginally stable tori, both I and £ are constant. Intending to modify the potential 
wn in such a way that its modified version w corresponds to the same family of 
equipotential surfaces as the general relativistic potential W, we have to require 

d r W _ C>rW 

d e W d e w' [ ' 

A natural modification of relation ([75]) is based on replacement of the Newtonian 
gravitational potential ifa by a new PN potential ip which fulfills condition ([77)) . 

In the SdS spacetime, the GR potential W is given by relation ([24|) . Because 
of the spherical symmetry of the spacetime, we assume ip ^ ip(0) and, applying the 
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Figure 13. Behaviour of the relativistic potential W in the equatorial plane of the 
SdS spacetime with y = 10 — 26 for two values of the specific angular momentum: 
I = £ mb = A.OAI (solid) and I = 1/2 (£ mb + £ ms , ) = 7.4 X 10 3 A/ (dashed), 
representing the most extended torus and a typical low density excretion torus. 



condition ([77)1 . we obtain an expression for the radial gradient of the PN potential ip 
in the form: 

w= 4-*^r (78) 

Now, it is clear that it is necessary to identify the Newtonian angular momentum 
per unit mass I with the GR specific angular momentum I = L/E to guarantee that 
the potential ip would be independent of the constant of motion, just as considered 
when defining ip by relation ([T]) in the general heuristic method. Then, of course, the 
integration of equation ([78]) yields tp in the form (fTTj) . i.e., identical with the heuristic 
method inspired by the test-particle motion. 



7. Conclusions 



It was shown that the PN gravitational potential ip defined for SdS spacetimes in close 
analogy with the well known Paczhski-Wiita potential reflects existence of the 
static radius, diverges at both the black-hole and cosmological horizons, and predicts 
locations of both the inner and outer marginally stable and marginally bound circular 
orbits at the same radii as those following from the full GR theory [7D]. The energy 
difference between the inner and outer marginally stable circular orbit, which plays a 
crucial role in the theory of thin discs, is very close to the relativistic result (for more 
details and comparison to the geodesic motion see [70]). Here, we have demonstrated 
that the PN potential ip can be well applied even in the case of thick discs orbiting 
SdS black holes. It provides exact determination of the equipressure (equipotential) 
surfaces governing the shape of toroidal discs in equilibrium configuration; to be 
precise, the shape of equipotential surfaces of the I — const tori are given by identical 
formulae when spherical coordinates are used in the PN approach and standard 
Schwarzschild coordinates in the GR approach. In a given SdS background, the 
marginally stable toroidal configuration with I — £ m b = const is the most extended 
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one. For y < 10~ 10 extension of marginally bound tori, expressed in terms of the 
gravitational radius r^ ~ M, falls with mass parameter M according to the law 
-Rout = r mb,o ~ A/ -2 / 3 . For £ = const < ^ mD , extension of the toroidal accreting 
configurations fastly falls even with extremely small decrease of t, similarly as the 
extension of excretion tori for extremely small increase of I above £ m b ■ 

We have shown that in both the PN and GR approaches there are very small 
differences when calculating potential barriers between the centre of the torus with 
£ = £mb — const and both its cusps. This indicates that physical properties of tori, 
predicted in both GR and PN approaches, could be very similar. 

We have tested physical properties of £ = const (marginally stable) toroidal discs 
in the case of adiabatic tori with fluid equation of state corresponding to the ideal 
gas. We have shown that the profiles of thermodynamic quantities (mass-density, 
temperature, pressure) are given by the profiles of cquipotcntial surfaces, and are 
almost identical in both the PN and GR approaches. Some differences appear in 
the central part of the tori, where temperature, mass-density and pressure given by 
the PN model are from a few percent to tens percent higher as compared with the 
GR model for current astrophysically relevant SdS spacetimes with y < 10~ 25 . The 
central temperature is given by the potential depth AW = W- m — W c only. On the 
other hand, the central mass-density (and pressure) is related the total mass of the 
disc. In the most extended discs with £ = £ m b, the PN and GR total masses of the 
disc are nearly equal, when a fixed central density is assumed, being a few percent 
higher in the GR model. Further, the non-relativistic fluid approximation can be 
conveniently applied for description of marginally bound adiabatic tori, since the total 
energy density exceeds the rest energy density in restricted central parts of the tori 
for a few percent only. 

Note that for y — > 0, extension of the £ = £„±, = const tori diverges and accretion 
tori with £ < ^ mD are relevant only. Precision of the PN model with £ — const < £ mD 
grows with £ — > £ mD . Therefore the results of the PN approach in Schwarzschild 
spacetimes [7T] can be applied for adiabatic tori. 

Extension of accretion (excretion) tori rapidly shrinks to the regions centered 
around the marginally stable circular geodesies r mSi i (r mS! o) when the torus parameter 
£ = const descends (rises) from the critical value £ — £ m b- 

In the moderate cases of accretion (excretion) tori, around currently 
astrophysically relevant SdS black holes, defined by the condition £ = const = 
l/2(£ ms .i+£ m b) (£ = const = l/2(£ m b+^ ms ,o)), with the fixed limiting mass m = M/10, 
significant differences appear. We find the central mass-density and pressure to be 
highly dependent on the black hole mass M and adiabatic index n for accretion 
tori; the differences of the central density and pressure obtained by the GR and PN 
models represent tens of percent, i.e., substantially more in comparison with the most 
extended tori. On the other hand, their central temperature is almost independent 
of the black-hole mass M. For excretion tori, the central mass-density is almost 
independent of the mass of supermassive black holes (M > 10 6 M©), while their central 
pressure and central temperature depend significantly on M, because of dependence 
of the potential depth on M in such excretion tori. The GR and PN models give 
quite negligible differences for these central values of thermodynamic quantities, not 
overcoming 10 -4 percent. This behaviour can be intuitively understood easily since 
the accretion tori under consideration are completely located in the strong gravity 
region near the black-hole horizon, while the excretion tori are in the region close to 
the static radius where the gravity is very weak and can be very well approximated 
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by the PN potential. The most extended £ = £ mh tori have only a small part located 
in the strong gravity region. 

We conclude that the PN modelling of tori in SdS spacetimes works very well 
in current astrophysically relevant situations with black holes of mass M < 10 10 M Q 
(corresponding to the cosmological parameter y < 10 -26 ) for tori with £ close to £ m b 
and for all excretion discs with £ m \, < £ < £ ms , , since physical properties of the fluid 
are properly reflected in these situations by the PN approximation. The PN approach 
can be useful even for the accretion tori in strong gravity regions, but the relativistic 
phenomena can shift the GR results for thermodynamic quantities by tens percent of 
PN results. 

We believe that the PN approach is useful in describing a wide range of 
astrophysical phenomena influenced by the observed relic cosmological constant. 
Especially, we can assume that the PN potential could be useful in developing models 
of excretion discs that could (potentially) serve as seeds of dark matter halos. 
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